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Abstract 

Dependences of low temperature behavior and anisotropy of various physical 
quantities for pure unconventional superconductors upon a particular form 
of momentum direction dependence for the superconducting order parameter 
(within the framework of the same symmetry type of superconducting pair- 
ing) are considered. A special attention is drawn to the possibility of different 
multiplicities of the nodes of the order parameter under their fixed positions 
on the Fermi surface, which are governed by symmetry. The problem of an 
unambiguous identification of a type of superconducting pairing on the basis 
of corresponding experimental results is discussed. Quasiparticle density of 
states at low energy for both homogeneous and mixed states, the low tem- 
perature dependences of the specific heat, penetration depth and thermal 
conductivity, the I-V curves of SS and NS tunnel junctions at low voltages 
are examined. A specific anisotropy of the boundary conditions for uncon- 
ventional superconducting order parameter near Tc for the case of specular 
reflection from the boundary is also investigated. 

PACS numbers: 74.25.Bt, 74.25.Fy, 74.70.Tx, 74.72.-h 
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I. INTRODUCTION 



The possibility of unconventional types of superconducting pairing in a number of heavy- 
fermion and high-Tc compounds, which follows from many experimental observations, results 
in much theoretical attention to a study of various physical properties of unconventional su- 
perconductors. If an unconventional superconducting order parameter A(p) vanishes at 
points or lines on the Fermi surface, low-energy quasiparticle excitations in the vicinities of 
these nodes may exist and arc known to give rise to power laws for low temperature depen- 
dences of the specific heat C(T), the thermal conductivity k{T) and some other quantities 
to be measured. Indices for these power laws depend not only on the dimension of nodes 
(that is whether there are points or lines of nodes) or, for example, on the strength of the 
impurity scattering. They differ also for the order parameters, which have zeroes of different 
orders at the nodes. Since always, even within the same symmetry type of pairing, there are 
basis functions which differ from each other by the multiphcities of nodes situated at the 
same lines and points on the Fermi surface, different indices may correspond to the same 
type of pairing. It seems to be quite important as different experimental consequences for 
superconductors with the same type of pairing become possible for different particular basis 
functions for the order parameter. For instance, let's consider (p^ — P^)-type of pairing in a 
tetragonal superconductor. For the basis function (pi — Py) and a spherical Fermi surface 
one gets at the Fermi surface A = Aq sin^ ^ cos 2(/7. Near fines of nodes (pi^2 — ±7r/4 and 
far from poles of the sphere this order parameter is proportional to {ip — ^pi,2) and near the 
intersections of lines at the poles one has A cc {(f — <fi^2)0'^- Other basis functions for the 
same representation of a symmetry group D^h and a spherical Fermi surface, for example, 
A = Aq sin^"* 6 cos^"^"*"^ 2cp (n, m = 0, 1, 2 . . .) have obviously different orders of zeroes at the 
nodes and, consequently, one obtains different low-temperature behaviors of various phys- 
ical quantities for these basis functions. The same arguments are applicable, of course, to 
the case of an arbitrary shape of the Fermi surface, for instance, to the cylindrical Fermi 
surface, when instead of the simplest representative for the order parameter A = Aocos2iy9 
for the case of (p^ — p^)-type of pairing, one can consider more complicated examples like 
A = Ao cos^""*"^ 2ip. Due to this fact an identification of the type of superconducting pairing 
(and frequently even a distinction between the contributions from lines and points of nodes) 
on the basis of experimental data for low temperature behavior of various quantities turns 
out to be ambiguous or at least more complicated as it is usually supposed. So, it is of 
importance to investigate the influence of the order parameter behavior near the nodes on 
some properties of unconventional superconductors. 

Having in mind this circumstance we calculate below the quasiparticle density of states 
at low energy and low temperature dependences of the specific heat ( for both homogeneous 
and mixed states ), the penetration depth, the thermal conductivity and the I-V curves for 
SS and NS tunnel junctions at low voltages, permitting the order parameter to have zeroes 
of different orders at the nodes on a Fermi surface. We consider also within this context an 
anisotropy of boundary conditions for an unconventional superconducting order parameter 
near for the case of specular reflection. Our results allow to carry out a qualitative and 
simple analysis of some experimental data, taking account of the possibilities mentioned 
above. We believe that such a consideration, which permits to discuss qualitatively experi- 
mental consequences for a wide set of the basis functions, is needed along with the detailed 
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quantitative (usually numerical) calculations of the effects for several currently preferable 
representatives from the set for the given compound. In particular, we discuss recent experi- 
mental data on the low temperature behavior of thermal conductivity for the heavy fermion 
superconductor UPts [00]. 



II. DENSITY OF STATES AND SPECIFIC HEAT 

If a superconducting order parameter doesn't vanish anywhere on the Fermi surface (as it 
is in the particular case of an isotropic s-wave superconductor), the density of states N{E) 
is equal to zero at low energies E < Amin and the electronic specific heat C{T) falls off 
exponentially with the decrease of temperature C oc e"^"*'"/-^, since only a small number 
of quasiparticles are activated above the gap at T <^ A^in- For an anisotropic pairing 
with nodes of the order parameter on the Fermi surface the low energy density of states is 
known to be nonzero. The contributions to N[E) at low energy come from quite narrow 
vicinities of the nodes, they dominates also in C{T) at low temperatures. Besides, in the 
case of unconventional superconductor the scattering on nonmagnetic impurities causes a 
substantial change in the low-energy density of states for various scattering strength and even 
at quite small impurity concentrations This underlies the low temperature behavior 

for various physical quantities including the specific heat. Below we consider only a pure 
unconventional superconductor which is relevant even to the case of resonance scattering 
at low impurity density, at least for the energies ujc ^ E <^ Tc above a critical energy 
cUc = (FAq)^''^ (where F is the normal state scattering rate in the unitarity limit). 

For a unitary gap matrix A(p) and a spherical Fermi surface the quasiparticle density 
of states in the case of pure superconductor is given by 

N{E) = 2y26{E-Ep)=NF f ^ , ^ , (1) 

P MiP)\<E An ^E^-\AipW 

where the integration is carried out over the momentum directions and Np is a density of 
states at the Fermi surface for normal metal. 

The electronic part of the specific heat may be written as 

E'^exp(E/T) rr^dNiE) , , , ^ 

' ' ' ^ ^ ^ ln(l + exp(-E/T)) + 



(exp(E/T) + 1)' dT 
E 



T(l + exp(E/r)) 



(2) 



One can easily estimate low energy and low temperature behavior of quantities N{E) 
and C(T) setting the particular form of the gap function only in the vicinities of nodes. For 
example, if the order parameter in the vicinities of lines of nodes takes the form (n > 0) 

|A(p)|=Ao|^-^or, (3) 
one easily gets the main contributions to N{E), C{T): 



3 



N{E) 
Np 



4h{n) {E/Ao 



l/n 



(4) 



C{T) = NF^w{n)T{T/A^ 



A/n 



Here the functions h{x), w{x) are given by 

h{x) 



r(l/2x) 



8v^a;r((x + l)/2)' 



(5) 



(6) 



wix) = 2h{x) (2 + ^) (1 - ^i^) r(2 + l/x)C(2 + 1/x), (7) 

where r(x), C,{x) are the Gamma function and the Riemann Zeta function respectively. 
Analogously, if near the lines of nodes one has {n > 0) 

|A(p)| = Aol^^-^ol", (8) 
the corresponding contributions to the density of states and specific heat are 



N{E) 
Np 



sin0o4vr/i(n) (E/Aq) 



l/n 



(9) 



C{T) = NpAiT sin eow{n)T (T/Aq)^/'' 



(10) 



There are obviously other basis functions belonging to the same representation of the 
symmetry group -D4/1. So, let the order parameter near the nodes on the Fermi sphere within 
the small angle regions \(f) — (pol < b, 9^ < a {a,h ^ 1) be 



mn2n 



\A{p)\ = AQ\^-^,re 

where m, n > 0, m + n > 0. Then one obtains from Eqs.(P, (^: 



N{E) 
Np 



m > n 



1—— 

(E/Aof'-f^Mm), 
{E/A.fl'' i (ln(Ao/i?) + c) h{n), m = n 



[{E/Aof''^h{n), 



n > m. 



(11) 



(12) 



CiT) 



NpTiT/Aof'^^Mm), 
NpT (T/Ao)'/" i (ln(Ao/r) + c) win) 



NpT{T/Aof''f^Mn), 



m > n 



m = n 



n > m. 



(13) 
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A role of logarithmic terms in Eqs.(|T2]), ([T3|) depends not only upon parameters E/ Aq, 
T/Aq but also on the magnitude of a constant c, which is formed by the contributions from 
the angular vicinity of the line of nodes, but not only very close to the pole. This constant 
may be calculated by making use of particular form of the basis function all over the Fermi 
sphere. For instance, in the case of {pi — p^)-type of pairing with A = Aq sin^ 9 cos 2{p one 
gets from (|l|) 

^^^^ -{A + \n{Ao/E))^, (14) 



where 



Nf 2Ao' 



1.38. (15) 



Here K{k) and E{k) are the complete elliptic integrals of first and second type. 

So, both terms in (|T^ may be of the same order and, as a rule, they both should be 



taken into account. Eq.(14) is in accordance with the second relation in (|T2D, since the 
order parameter A near the nodes 6o = 0, ipi^2 = ±7r/4 takes the same form as in (|TT|) with 
n = m = 1, and the logarithmic term in ([T^ ) contains respective contributions from all four 
nodes (that is from two lines of nodes at two poles of a Fermi sphere ). As in this example, 
the behavior of the order parameter of the form (|Tl]) usually corresponds to the intersection 
of two lines of nodes at the pole (each being of the form |A(p)| = AqIv? — v^ol"^ near the 
line and far from the pole). One can easily seen that under the condition n > m the part 
coming from the point of intersection dominates or is at least of the same order as the total 
contribution to the quantities. Note, that in the particular case m = Eqs.(0), (|T3p at 
n > m correspond to the contributions from the point node which is described by Eq.(P) for 
^0 = and the substitution n — > 2n. 

At last, for the order parameter which may be represented near the point node as follows 

\A{p)\ = \Aoi{ip - ipor + Ao2{9 - 9o)% m,n>0, m + n>0, (16) 



one gets 

N{E) 
Np 



ocsin^o^'/"^'/", C(T) ocsin^or'+i/"+'/™. (17) 



Summing up the contributions from all vicinities of the nodes of the order parameter, 
one finds total low energy density of states and low-temperature specific heat. 

In discussing lines of nodes of the order parameter one usually considers zeroes of the first 
order at the nodes which correspond to relation (|^) or (|]) with the particular value n = 1. 
Then, according to the Eqs.(^, (p!0|) one gets a quadratic low-temperature behavior for the 
heat capacity C{T) oc T^. Analogously, for the points of nodes it is frequently supposed 
n = m = 1 in (|T6|) and then one obtains from (0) the well-known result C{T) oc T^. As it is 
seen from the more general relations (H), ([To|) , (0), the possibility for distinguishing between 
the contributions from points and lines of nodes becomes substantially more restricted, and 
indices of the power laws are allowed to be fractional, even for integer values of n, m. For 
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instance, for an experimental determination whether there is a contribution from a cubic 
point of nodes (m = n = 3, when corresponding term in the specific heat is Ccp oc T^/^), or 
from a first order hne of nodes (when one has Cfi (x T^), quite a high accuracy is needed 
as compared to the distinguishing between the heat capacity dependences Cfp oc and 
Cfi cc T^. Note that a contribution to the specific heat from zeroes of second order at points 
(quadratic points) is proportional to precisely as from the first order line of nodes. In 
the case of fractional values of n, m as well as for the dependences like |6' — the question 
about the physical origin for the nonanalitical behavior of the order parameter near the 
nodes appears. 



III. LOW ENERGY DENSITY OF STATES AND LOW-TEMPERATURE 
SPECIFIC HEAT FOR THE MIXED STATE 



One of the most important qualitative features of an unconventional superconductor 
in a mixed state is that due to the existence of a supercurrent induced by a magnetic 
field, the quasiparticle states with negative energy (which are obviously occupied) appear 
for momentum directions within the vicinity of nodes of the order parameter, as it follows 



from the relation E{p) = \j i{pY + |A(p)P + pvg 0. The density of states for such a 
superconductor changes substantially at sufficiently low energies. In particular, the density 
of states at the Fermi surface differs from zero under these conditions and after spatial 
averaging over the mixed state is proportional in the simplest case to \/B, resulting in the 
characteristic magnetic field dependent term in the specific heat C cx T^/B 0. Recently the 
term of such kind, which is specific for unconventional superconductors, has been observed 
experimentally for YBCO From these measurements a quite strong upper limit on 

the minimum gap value follows if an anisotropic s-wave superconductivity is assumed. We 
examine below the infiuence of the particular forms of basis functions, which belong to 
the same symmetry type of pairing, on the magnetic field dependence of the density of 
states and heat capacity (under the conditions similar to those considered in |]^). It will 
be shown that an accurate experimental determination of the index value in this power law 
dependence permits to describe the behavior of the order parameter in the vicinity of nodes 
more carefully. The spatial distribution of the density of states around the vortex is also 
considered. 

Let a strong type two pure superconductor with a cylindrical Fermi surface be at low 
temperature in a mixed state under the applied magnetic field, which satisfies the condition 
Hci B <^ Hc2 and is directed along the principal cylindrical axis. For this inhomogeneous 
state the local density of states must obviously manifest a spatial dependence and, generally 
speaking, is formed by both a superconducting region outside the vortex cores and the 
quasiparticle states localized inside the cores. Since the total volume occupied by the vortex 
cores is proportional to the magnetic field, the contribution to the density of states from 
the states localized within the vortex cores is a linear function of B and the corresponding 
contribution to the specific heat cx TB/Hc2 [llO|. The quasiparticles at distances r ^ 



from the vortex core may be considered quasiclassically as having the energy 



E{p,r) = J^py + \A{p)\^+pv,{r), (18) 
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which depends upon the distance from the vortex core along with the superfluid velocity 
(in the simplest case of circular supercurrents one has Vg = e^Ki{r / X) /2m(.\) . Then the 
corresponding density of states may be written as follows 



N{E,r) = 2j ^^S^^e{p) + \MPW + Pvs{r)-Ey (19) 

From the estimate \pVs\ ^ pp/mer <^ vp/^o ~ Aq which is valid at distances r ^ 
it follows, that a nonzero contribution to the density of states for the zero quasiparticle 
energy appears only for unconventional superconducting pairing for momentum directions 
within the narrow vicinities of nodes of the order parameter A(p) ( Aq denotes throughout 
this paper the maximum value of the order parameter for the homogeneous superconducting 
state at low temperature and for any type of pairing). Due to this fact in considering the 
density of states at quite low energy one can describe the contributions from each line or 
point of nodes separately. 

The location of a line of nodes of the order parameter, which is oriented parallel to the 
principal axis of the cylindrical Fermi surface, may be described by a constant polar angle ipi 
of the cylindrical coordinate system, corresponding to the momentum direction Pf,±(0 to 
line. In integrating in (|T9|) one can put with a good accuracy pVg = PF,L{i)vs due to a quite 
small parameter Vg/vp- Suppose that the quasiparticle energy ^(p) for the normal metal and 
the order parameter A(p) near nodes depend upon the different momentum components - 
the magnitude of the momentum component Pr and the direction of the momentum Pp_i 
(which is perpendicular to the principal cylindrical axis) correspondingly. Then one obtains 
after integration over the energy ^{pr) {dC, = vpdpr) the following contribution from the line 
of nodes 

Here G(x) is the step-like function, Np is the density of states at the cylindrical Fermi 
surface for the normal metal and Qi is the narrow angular region in the vicinity of the line 
of nodes, where the expression under the square root sign in (pOf) is positive at quite low 
energies in question. 

If in the close vicinity of the nodes the order parameter takes the form 

|A(<^)| = Ao|^-<^d", n>0, (21) 
one gets from (^Oj) after the integration over 



NiE, r) ^ ^^^Ht f ^-"^f"-^"^ ) e(^; - pMr)). (22) 



nx/vFr - + — 
2n 



According to this relation, the narrow angular region in the momentum space situated 
near the line of nodes on the Fermi surface, contributes to the local quasiclassical density of 
states at sufficiently low energy for quite a wide region of orientations in the coordinate space. 
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For instance, at zero energy the density of states formed by the hne of nodes is nonzero for 
the whole spatial region defined by the condition Pf,±(0'^s('") < 0) though it naturally falls 
off with the increasing distance from the vortex core (along with the supecurrents induced 
by the magnetic field). So, in the presence of several lines of nodes the contributions are 
spatially superimposed and due to this fact the locations of a maximum value of the total 
density of states in the coordinate space (at the given distance r from the vortex core) may 
substantially differ from the locations of corresponding maxima of the partial contributions 
from each separate line of nodes. 

In the case of {p^ — Py)-pairing for the tetragonal superconductor with cylindrical Fermi 
surface, one can consider four similar lines of nodes for the superconducting order parameter, 
which are situated at the angles ipi = 7r/4+(/ — l)7r/2 (/ = 1, 2, 3, 4). Let the angle specifies 
the relative orientation of the vector r and the crystalline axis x, and Vs{r) = Vs{r)e^. Then 
the sum of all four terms of the form ( ^2]) results in the following spatial dependent density 
of states 

iV(E,r,0) = 

j ^(E + sin(0 + y;O) "e(E + sin(0 + (/;;)), (23) 



2nJ VFKi{r/\) 



^ /I 1 \ I 2AAo 

where the dimensionless quantity E = 2XE /vpKilr/ X) is introduced. 

Since at large distanses from the vortex core the anisotropy of is small, we con- 
sider further only the distances r ^ A. At the zero energy one gets from here A^ cx 
Npi^o/rY^'' (l cos(0 + 7r/4)|i/" + | sin(0 + 7r/4)|i/") and in the particular case n = 1 the 
spatial angular dependence of A^ is simply A^ oc |cos0| for —tt/A < (p + -Kin < it /A and 
A^ oc I sin 01 for n/A < (p + Trm < 3tt/A (m = 0, 1). Hence, the maximum value of the density 
of states at given r lies at = 0,±7r/2,7r and minimum - at = ±7r/4, ±37r/4. At the 
same time the maximum of the separated contribution from each line of nodes lies in the 
coordinate space at the direction, where the superfluid velocity is opposed to the momentum 
direction to the corresponding line on the Fermi surface (i.e. all four these maxima lie in 
the directions = ±71/4, ±3-71/4). Furthermore, in the particular case n = 1 and for E > 1 
(e.g. at the distances r > (Tc/ E)^q) all spatial dependence of the density of states (both 
from the distance and the direction) disappears and one gets from ( PBD A^ ~ {E / Aq)Np. 
It follows from here that the greater the energy the smaller the spatial region where the 
anisotropy and inhomogeneity of the density of states manifest themselves. So, for ~ 
one has the anisotropy only within the region r < C,o, where the local description based on 



the relation (|I8D for the quasiparticle energy, can't be applied and more general approach is 
needed [jl3[ . Therefore, within the framework of the approach used above one may consider 
only the energies E -C Tg. It is worth noting that in contrast to the particular case n = 1, for 
n ^ 1 the angular and distance dependences of the density of states at large distances don't 
disappear completely though become quite weak. The spatial dependence of the density of 
states (0), including its fourfold symmetry, is accessible for studying to scanning tunneling 
microscopy experiments, which now have a high spatial resolution as a valuable complement 
to the high energy resolution of tunneling spectroscopy (see, for example, |]TT],|T2[). 
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As the density of states at zero energy is finite, then one gets from (0) the conventional 
(hnear in T) expression for the low-temperature specific heat C{T) = tt'^TN{0)/3. Due to 
the inhomogeneity of the state the spatial average of the density of states iV(0) is presented 
here. For the magnetic fields within the interval Hd ^ -B <C Hc2 one substitutes A/r instead 
of Ki{r/X) in (0), (^) and then the spatial averaging over the region C,o <ti r <ti R (where 
R ~ ^o{Hc2/ BY^'^ is the distance between vortices) gives rise to the following estimate 

/ VF V'"" f B 



Here depends only upon the index n. Note, that in forming iV(0) just the distances 
r ^ ^0 are found to be important for n> 1/2. For these values of n the contribution to the 
specific heat oc T{B / Hc2Y^'^^ is more essential, than the linear magnetic field dependence, 
coming from the interiors of the vortex cores. 

According to (^), the index l/2n of the power law for the magnetic field dependence of 
the low-temperature specific heat for the unconventional superconductor in the mixed state 
characterizes not only the existence of nodes of the order parameter on the Fermi surface 
themselves, but also the behavior of the order parameter in the vicinities of the nodes, for 
example, of a form ([21|). 

IV. PENETRATION DEPTH 

It is well known that deviation of the penetration depth at low temperatures from its 
zero temperature value is exponentially small in the case of a superconductor with finite gap 
and manifests a power law temperature dependence for a superconductor with nodes of the 
order parameter on the Fermi surface p!5|-p!9|]. We are interested here in the indices of those 
power laws for temperatures (FAq)^/^ < T ^ Aq in the case of pure superconductors with 
unitary order parameters having various multiplicities of the nodes. As a starting point one 
can use the following expression for the eigenvalues of the penetration depth tensor: 

^ W /• d'S „ |A(p)|'= 



where index i denotes the direction of the supercurrent and vp^i is the i-component of the 
quasiparticle velocity at the Fermi surface. Applying the Poisson formula to the sum in 
(E3), we get for the deviation Aj(T) from Aj(0) at T ^ Tc : 



J Vp 1 r, 



5A,(T) J VF 



m=l,2,... 




Ai(0) fd'S 



(26) 



Vf 



where Ki{x) is the modified Bessel function. 

For T -C Aq the main contribution to the integral over the Fermi surface in the numerator 
comes from narrow regions near nodes of A(p). As in the previous section we consider a 
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tetragonal superconductor with cylindrical Fermi surface and the order parameter having 
four lines of nodes. Then it follows from ( p6l) for the in-plane penetration depth 



A||(0) Jo 2,r„ij3j '\ T J T 

where the integration is carried out over the angle in the basal plane. 

Assuming relation ( pT]) for the order parameter in the vicinities of lines of nodes we get 
for the leading correction to the penetration depth at low temperatures 

C^inUir, (28) 



A||(0) ' ' \A, 



^0 



where 



^ ^ Tin \ 2n J \2n) w}''^ 



In the particular case n = 1 one obtains from here 5A|j(T)/A||(0) = 21n2 (T/Aq). 

V. THERMAL CONDUCTIVITY 

The investigation of the thermal conductivity in a superconducting state is an important 
experimental probe of the order parameter structure. 

Thermal conductivity for an isotropic s-wave superconductor is known to fall off expo- 
nentially at low temperatures due to the presence of a finite superconducting gap for all 
directions on the Fermi surface. For unconventional superconductors a temperature depen- 
dence of the electron thermal conductivity at low temperatures T ^ Aq is described by 
the power laws and proves to be associated with several factors. Besides the order param- 
eter structure, it is substantially influenced by the strength of impurity scattering. For the 
thermal conduction components directed to the nodes of the order parameter one usually 
obtains in the case of weak scattering processes (Born approximation) a linear temperature 
dependence, as in the normal state. Other components vanish with the temperature accord- 
ing to the higher power laws. In the case of resonance impurity scattering the situation is 
more complicated and the region of low temperatures T -C Aq is naturally devided into two 
parts. For temperatures T < cj^ = (r^o)^'^^ the existence of the quasiparticle states bound 
to impurities is of importance. The correspondent thermal conductivity at those temper- 
atures has also linear temperature dependence albeit with reduced density of states. By 
contrast, for the temperatures Uc ^ T <^ Aq one can disregard the bound states and then 
one obtains, generally speaking, the higher power laws for all the components of thermal 
conductivity. Since experiments have shown that in several heavy fermion superconductors 
the thermal conductivity vanishes with higher powers of T at low temperatures (see, for 
example |^) the theory for thermal conductivity in unconventional superconductors has 
been developed and examined in detail pll-[2^j^,p5|-|3l| . The experimental results for UPt^ 
are in correspondence with the strong scattering from the impurities, close to the unitarity 
limit. Some physical discussion of the strength of the impurity scattering in UPt^ is in 
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Recent experimental results for UPt^ obtained in |l|J^ have drawn much attention 
[l], as they gave some further useful information about the order parameter structure 
in this heavy fermion hexagonal superconductor. In particular, it has been observed that 
the anisotropy ratio k^/ Hb of the thermal conductivity components in UPts does not vanish 
at low temperatures (down to Tc/10). The temperatures, at which one can disregard the 
influence of impurity bound states on the transport coefficients, depend upon the impurity 
concentration. It has been stressed earlier that for superconducting UPt^ this influence 



usually may be neglected just down to Tc/10. It is conflrmed by the absence of a linear term 
in the temperature dependence of the thermal conductivity observed in for temperatures 
down to Tc/10 as well. 

Based on this fact below we give a simple qualitative examination of the thermal con- 
ductivity at low temperatures for pure unconventional superconductors at temperatures 
cUc ^ T <^ Aq, separating the contributions from points and lines of nodes (and from the 
rest of the Fermi surface, when it is needed) and permitting the order parameter to have 
zeroes of different orders at the nodes. Such a consideration results in a simple condition for 
the realization of a flnite anisotropy ratio Hc/ i^b at low temperatures. We believe that our 
simple consideration of the thermal conductivity at low temperatures tUc ^ T <^ Aq, which 
permits to discuss qualitatively a wide number of basis functions for various types of pairing, 
complements in some aspects several recent quantitative investigations of this problem (see, 
for example, ||28|j30|j3l|] ) . On the basis of self-consistent description of the impurity scattering 



in unconventional superconductors, the numerical calculations are performed there for all 
temperatures and for several particular basis functions, which are considered now to be of 
special interest for the study of superconductivity in UPt^ and in HTSC as well. 

So, lets neglect the impurity bound states and consider only essentially anisotropic super- 
conducting states which meet the condition J2p ^{p) / {E^ ~ E^) = 0. Then one has for the 
thermal conductivity of unconventional superconductor in the case of resonance scattering 
(see e.g. mm) 



kn{T,) tt'tJo VT 



dn%E) 
dE 



giE)\' r dnp,p,^/E'-\Aip 
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Reg{E) J\A{p)\<E Att p^ E 



(29) 



Here ki^{Tc) is the thermal conductivity in the normal state at T = Tc, n^{E) is the equilib- 
rium Fermi distribution function for excitations, and the function g{E) is deflned as follows 

g{E) = / , 30 

^ Je2- |A(p)|2 



One can easily see that due to the factor dinP^E) / dE at low temperatures T <^ Aq 
only low energy behavior of the integrand in (|29| ) is of interest (for E ^ T). Then only 
narrow regions near the nodes of the order parameter are essential for the integral over the 
momentum directions in (^) and for the real part of the function g{E). But, generally 
speaking, it is not the case for the imaginary part of g{E) at low energy (see (|30[) ). Due 
to this fact and contrary to the case of speciflc heat, for the thermal conductivity at low 
temperatures T <^ Aq contributions from all the Fermi surface - not only from the narrow 
vicinities of nodes of the order parameter on the Fermi surface - may be of importance. 
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Nevertheless the contributions can be divorced from each other and specificated, and in the 
case of zeroes of higher orders at the nodes only those narrow regions become dominating. 
Note that regions far from nodes, as it is seen from (pO]), result in the imaginary part of 
g{E) at low energy which is simply a linear function of energy. 

Let's consider, for instance, the order parameter, which behaves near the line of nodes as 
\A{p)\ = Ao\e - 7r/2|" {\e - tt/2\ < 1, n > 0) and near the point node as |A(p)| = Aq^"" 
< 1, m > 0). Then one gets at low energy Re g{E), Im g{E) oc (E/Aq)^/", (E/Ao)^/™ 
respectively, but in the case of imaginary part of g{E) this estimate is valid only if n > 1, 
m > 2. For values m < 2 and n < 1 one obtains for the main term Img{E) (x E which 
results from regions lying far from the nodes of the order parameter on the Fermi surface. 
So, both Reg{E) and Img{E) should be taken into account in considering the thermal 
conductivity at low temperatures. In the case n > 1, m > 2 only the behavior of the order 
parameter near nodes turns out to be essential; for n < 1, m < 2 contributions from all the 
Fermi surface are of importance. For specific and quite important particular values n = 1 
and m = 2 one gets apart from the linear terms in E for Reg{E), Img{E) also a contribution 
from the vicinity of nodes Img{E) oc {E / Aq)\h{E / Aq). As in the case of specific heat a 
pure logarithmic approximation is a very strong restriction on the magnitude of energy and 
usually logarithmic factors must be taken into consideration together with constants. The 
anisotropy of the thermal conduction is defined by the integral over momentum directions 
in (|29|) as well. As a result one gets for temperatures oJc^T <^ Aq: 

K,, = T5(ln2(T/Ao) + hi ln(T/Ao) + /.a), n = 1 + 
[ T2/"+3, n < 1 

r m>2 

+ pA T3(ln2(r/Ao)+j9,iln(T/Ao)+p.2), m = 2 (31) 
[ T^, m< 2, 



yy — 



rpl+2/n^ n > 1 

T^\n\T/Ao) + UHT/ Ao) + h^), n = l + 



PA T4(ln^(r/Ao)+p.iln(T/Ao; 

rp3+2/m 



m> 2 
+ Px2), m = 2 
m<2. 



(32) 



The coefficients L, In, li2 and P, pn, pi2 {i = x, z) correspond here to the contributions 
from the line of nodes and the point node respectively. 

It follows from Eqs. (|3lD , (p2|), that under the condition m = 2n the main contribution 
to K,zz at low temperatures originates from the point node, while in the case of k^x the 
contribution from the line of nodes dominates. Taking this fact into account one obtains in 
qualitative agreement with experimental results the finite (temperature independent) value 
for the anisotropy ratio k^z/ i^xx at low temperatures under the condition m = 2n. It is worth 
noting that in the case of hexagonal superconductor an important particular example for 
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the type of pairing with m = 2n (and integer n, m) is a £'2?i-representation (in the simplest 
case m = 2n = 2). At the same time in the simplest case of i^ig-representation one should 
put n = m = 1 and then it follows from Eqs.(pl]),(p^, that the ratio Kzz/k,xx manifests only 
feebly marked logarithmic temperature dependence at low temperatures, which possibly can 
not be determined now experimentally within the temperature interval uJc ^ T <^ Aq. 

From this qualitative discussion one can arrive at the conclusion, that the both repre- 
sentations Eig and with the simplest particular basis functions do not contradict to the 
finite value of the thermal conductivity anisotropy ratio at temperatures ^ T -C Aq. In 
more complicated cases the anisotropy ratio k^z/ i^xx in fact doesn't depend upon tempera- 
ture within the temperature region considered, if the indices n and m satisfy the condition 
m = 2n. The presence of both line of nodes and point nodes is essential in order to form 
the observed temperature dependence of the thermal conductivity. 



VI. QUASIPARTICLE TUNNELING AT LOW TEMPERATURE AND VOLTAGE 

For isotropic s-wave superconductors a small number of quasiparticles activated above 
the gap at low temperatures T ^ Aq gives rise to the junction conductance which falls 
off exponentially with decreasing temperature. At zero temperature and for the externally 
applied voltage V < Ao/e for NS junctions (and V < (Ao+ + Ao-)/e for SS junctions) the 
quasiparticle current across the tunnel junctions takes place only for superconductors with 
anisotropic order parameter. In the case of very small value of voltage {V <^ Aq) the current 
occurs only for the superconductors with nodes of the order parameter on the Fermi surface. 
Below we show that respective I-V curves for these junctions essentially depend on the 
behavior of the order parameter in the vicinity of nodes, in particular, on the multiplicities 
of the nodes. 

If the permanent voltage V is applied to the tunnel junction between two metals, then 
one gets the following expression for a dissipative contribution to the current across the 
junction in the lowest order in the junction transparency D{p) [p = Pf/\Pf\', see e. g. 



my- 



JN 



Here g±{E,p^) are the normalized densities of states (for fixed both energy and the mo- 
mentum direction) of two metals in the vicinity of a tunnel barrier, the index -|-(— ) labels 
the right (left) half space with respect to the boundary plane, and Vx is the Fermi veloc- 
ity component along the normal to the plane interface n \\ Ox. The integration in (|33D 
is carried out over the part of the Fermi surface with > 0. The relation between the 
incident and transmitted Fermi momenta ( that is between p_ and p+ ) is as follows. The 
components parallel to the specular plane interface are equal to each other, while the values 
of their normal components are determined by the Fermi surfaces of corresponding metals. 
Naturally, in the particular case of identical superconductors with the spherical Fermi sur- 
faces the total incident and transmitted momenta are equal to each other p_ = For an 
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anisotropically paired superconductor the tunneling density of states is essentially depen- 
dent upon the orientation of crystalline axes of the metal relative to the interface. Below we 
consider, for simplicity, only those particular crystal orientations, when the order parameter 
is not suppressed at the specularly reflecting interface as compared to its value in the depth 
of pure superconductor. These orientations are determined by the equation A(p) = A(p), 
where p [p) denotes the direction of the incident (reflected) electron momentum (see e.g. 
Then the following expression for the density of states is valid 



\E\Q{\E\ - \A^{p^)\) 



- \A 



±[P±} 



Under the condition T ^ eV Eq. (|33|) is transformed after substitution (^J) into 
the expression 



(34) 



to 



= 2eV / / 

Jvj:>o Jo 



duo ^ 



A-{p_)/eV 



A+{p+)/eV 



A_(p_)|V(ey)V(l - uof - |A+(p+)|7(er)^ 



X D{p_)v,^{p_ 



(35) 



One can easily see that for SS junction in the case \eV\ <^ Ao± the contribution to the 
integral in (|35D comes entirely from the narrow vicinities of the directions to the common 
nodes of two order parameters, as if they were drawn on the same Fermi surface in taking 
account of the relative orientation of the superconductors. These directions are defined by 
the relation |A_(p„)| = |A+(p_j^)| = 0. Let the latter condition results in the momentum 
direction p^ which corresponds to a point of the intersection on the Fermi surface of two lines 
of nodes of the order parameters A_ and A+ respectively. We suppose that the behavior of 
the order parameters close to this point may be represented as follows 



|A-(P) 



A, 



o-|7- 



|A+(p)| = Ao+17+1 



n, m > 0, 



(36) 



where 7^: are the angles which are counted from p^ in the directions perpendicular to the 
lines of nodes of the corresponding order parameters A^. Substituting (|36|) into ( ^5]) we 
obtain in the particular case of spherical Fermi surface the leading order contribution to the 
quasiparticle current for eV ^ Aq : 



3n 



a- 



Vp\ smxi 



MPi)D{pMeV)\eV/Ao\ 



(37) 



Here a is a numerical factor and x is the angle between two lines of nodes in the point of 
their intersection. The angle x is supposed to obey the condition x ^ {eV/ AqY^"^ ( for 
m > n). In the opposite limit the lines coincide and then one gets Jn oc V^^+(^/"). Summing 
up the contributions from all intersection points one finds 



JN 



R 



-V\eV/Ao 



■N 



(38) 
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where Rn is the normal-state resistance of a tunnel junction and the numerical factor a 
depends on the relative orientation of superconductors. 

We see that the index of the power law ( |38D for the I-V curve in the case of SS junction 
at law voltage and temperature depends essentially upon the multiplicities of nodes of the 
order parameters. In the particular case n = m = 1 the result ( p8| ) is reduced to that 
obtained earlier in . 

The quasiparticle current for the case of NS junction may be obtained by substituting 
A_|_ = into {EE) and integrating over u 



jN = 2eV / . 1 



|A(p) 

{evy 



-0 1 



A(p) 



eV 



(fS 



{2'nfvF' 



(39) 



If the externally applied voltage is small enough, so that \eV\ -C Aq, the contribution to jtv 
in (^) comes from narrow vicinity of nodes of the order parameter. One can easily verify 
that if the order parameter near the lines of nodes on the spherical Fermi surface has the 
form |A(p)| = Ao|6' — 6'o|" or |A(p)| = AqIv? — v^oT; it follows from (^) for the quasiparticle 
current across NS junction 



ll/n 



Analogously, in the case of the order parameter behavior |A(p) 
Ao2(^ — 6qY\ close to the point node in the direction Pq, one gets 

V , , ,1 ,i 
jN OC — — |eV^/Ao|"^-. 

At last, for the order parameter behavior near the pole |A(p)| : 
contribution to the quasiparticle current is given by 



(40) 

= |Aoi(v9 - + 
(41) 

Ao|0 - 00 r^^" the 





V 


eV 
Ao 




V 


eV 
Ao 




V 


eV 
Ao 



l/ni 



l/m 



l/n 



In 



Ar 



eV 



+ c 



m > n, 



171 = 71, 



m < n. 



(42) 



The results ( p]| ) and (^) in the particular cases n = 1 and m 
those obtained earlier in ||33ll. 



n 



1 are reduced to 



VII. ANISOTROPY OF THE BOUNDARY CONDITIONS FOR A 
SUPERCONDUCTING ORDER PARAMETER 

Many aspects of the problem of boundary conditions for the order parameter of 
anisotropically paired superconductors have already been studied in the literature (see, 
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e.g. |]3^-p^, ^ , p7[] and references therein). Below we discuss one characteristic feature of 
the boundary conditions near Tf. , which has been recently noticed in for the case of 
specularly reflecting boundary. For simplicity, let's consider a singlet pairing and a one com- 
ponent anisotropic order parameter in the vicinity of the plane specular and fully reflecting 
boundary at x = 0. For any crystal orientation with respect to the interface the boundary 
condition may be written in the form qri'{0) = ri{0). It turns out that with the aid of general 
symmetry arguments and, for example, irrespective of the shape of the Fermi surface of the 
given symmetry, one can unambiguously fix only those crystal orientations for which the 
parameter q is equal to zero and infinity. The value q = oo corresponds to the boundary 
condition ^^'(O) = coinciding with the condition for an isotropic s-wave superconductor at 
an impenetrable boundary, where the relation rj^O) = rj^ is valid. Here r^oo is the value of 
the superconducting order parameter far from the boundary. For an anisotropically paired 
superconductor one has q = oo only for those orientations of the boundary for which the 
relation A(p) = A(p) holds. The value g = is realized under the condition A{p) = —A{p), 
in which case the boundary condition takes the form 77(0) = and the order parameter is 
fully suppressed at the interface. As earlier p (p) denotes here the direction of the incident 
(reflected) electron momentum. For the point symmetry group D^h of the crystal there 
are five planes of symmetry, and the order parameter transforms according to one of the 
irreducible representations of this group. So, if the normal to the boundary is perpendicular 
to the symmetry plane and the character of the irreducible representation is 1 (-1), the 
boundary condition is V(0) = (^7(0) = 0)- 

For other orientations the detailed form of g as a function of crystal orientation relative to 
the boundary plane for the given type of pairing essentially depends on the particular form 
of the basis function A(p) in the depth of superconductor. Since the characteristic value for 
the parameter q is ^(T) one can divide all orientations into two regions where q < ^(T) and 
q > C,{T) respectively. At first sight the characteristic scales for these two regions could be 
comparable on the order of values. However, it is not the case near Tc and may be valid only 
for low temperatures. E.g. for the particular basis function pl^ — ( pertaining obviously 
to the pairing symmetry of the type p^Q — p^q) and the spherical Fermi surface, the values 
ri{0) and r/oo are found to be of the same order of magnitude (and hence ^(T) < q) only 
for particular orientations of the normal to the boundary n within narrow angular intervals 
A0O ~ A^o ~ {Co/CiT)Y^'^ around the crystal axes xo,yo,zo. For other crystal orientations 
one has 77(0) ~ r]oo^o/^(T) -C r/oo- Thus, near Tc the d-wave superconducting order parameter 
proves to be strongly suppressed in the vicinity of the insulating barrier, except for the small 
part of orientations within the narrow angular intervals mentioned above. From this point 
of view the qualitative difference between the boundary conditions for specular and diffusive 
interfaces near Tc manifests only within those narrow angular intervals. Indeed, in the case 
of diffusive scattering at the interface one has q ^ ^0 irrespective of the orientations and, 
hence, near Tc the parameter 77(0) ~ ^oVoo/^iT) <^ t]^ for all crystal orientations. So, near 
Tc narrow peaks ^0 ^(T) < g of the parameter q, which occur in the case of specularly 
reflected boundaries, are simply cut if the boundary is diffusive. 

It is of interest to consider whether the choice of a different basis function pertaining to 
the same symmetry type of pairing can result near Tc in the expansion of the angular regions 
where the parameter q has large values ^(T) < q. Below we show that it is not possible for 
quite a wide set of the bases functions. In particular, for the bases functions with higher 
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multiplicities of the nodes on the Fermi surface the angular regions with ^(T) ^ q prove to 
be more narrow, not wider. 

Microscopic consideration based, in particular, on some variational procedure, leads to 
the following expression for the parameter q near [^ , ^ ,^ : 

y (43) 



X 

Here integration over the Fermi surface is confined to its part with f ^ > and the following 
notations are introduced: Vx = vfx/vf, = vp/Tc- The superconducting order parameter 
for the inhomogeneous state near is supposed to be represented in the conventional form 
A(p,r) = ip{p)r]{r). 

For isotropic superconductors the basis function ip doesn't depend upon the momentum 
at all and one gets from (|i3|) the equality q = oo. For anisotropically paired superconductors 
there rule, only several isolated and governed by the symmetry crystal orientations 

relative to the boundary for which q = oo and, hence, ip{p) = ip{p) for all p. The crystal 
orientation with respect to the boundary may be characterized by the unit vector n which 
describes the orientation of the normal to the boundary relative to the crystalline axes 
xo.uo.zo. Denoting by no the orientation corresponding to the value g = cxd, we are interested 
now in the characteristic scale of the angular regions around the direction no where ^ (T) ^ g, 
i.e. the parameter q has quite a large value yet. 

It follows from Eq . (|43D , that the relation ^(T) ^ g is valid under the condition 

Since in the case n = no one gets ip{p) = ip{p) (in the particular case when the boundary 
coincides with the symmetry plane of the Fermi surface, one has p = p — 2(pn)n), in 
the narrow vicinity of this direction one obtains ip{p) ~ ip{p) + (<5^Vn)'?A(p)|n=no- Here 
3n = n—riQ and one has (n, no) ~ |(5n| for the angle between vectors n and no. Substituting 
these relations into (^4]), we obtain the following restriction for the angle (n, no): 



^^^^^ (/..>0 \{eSr^V^)^{p)W=.^^vJ^S) 

where e^n is the unit vector along the direction 5n. 

We see that the appearance of the small parameter (^o/^(^))^^^ is a quite general char- 
acteristic feature for the estimation of the angle (n, no), which is independent upon the 
particular form of the basis functions. In order to estimate the other factor by which the 
first one is multiplied in (|45|) , let's consider the following set of the basis functions for the 
tetragonal unconventional superconductors 

V'(p) = Aotanh{/i(p^Q-i}y"} , n=l,2,... , (46) 
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where /i and n are the parameters. All functions with odd values of n belong evidently 
to the same type of pairing as the function (p^Q — p^o)- Analogously, the functions with 
even values of n pertain to the unit representation of the point tetragonal group D^h (then 
respective nodes of the order parameter, of course, are accidental and not associated with 
the symmetry type of pairing ). 

The substitution of (^61) into 
intervals around no, where ^{T) ^ q: 

where v = max(l, fi). 

It follows from here that the angular intervals can't be essentially more than {^o/^(T)y^'^y 
but may be more narrow, for instance, in the case of higher multiplicities of the nodes of 
the order parameter. 

VIII. CONCLUSIONS 

Density of states at low energy, the specific heat in homogeneous and mixed states at 
low temperatures, the low temperature behavior of the thermal conductivity, the penetration 
depth and I-V curves for the quasiparticle current at low voltage have been examined for 
pure anisotropically paired superconductors having various multiplicities of the nodes of the 
order parameter. A specific anisotropy of the boundary conditions for unconventional su- 
perconducting order parameter near for the case of specular reflection from the boundary 
was also investigated. Since the multiplicites may be different within the same symmetry 
type of pairing, the problem of an unambiguous identification of a type of superconducting 
pairing on the basis of corresponding experimental results becomes more complicated and 
should be considered taking into account the results obtained above. 
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